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ИСПОЛЬЗОВАНИЕ ПРОГРАММНЫХ СРЕДСТВ ПОДДЕРЖКИ 
ПЛАНИРОВАНИЯ ЭКСПЕРИМЕНТА ДЛЯ ОПТИМИЗАЦИИ ПАРАМЕТРОВ 
ГЕНЕТИЧЕСКОГО АЛГОРИТМА АППРОКСИМАЦИИ 


Описывается экспериментальная оптимизация скорости сходимости генетического алгоритма аппроксима- 
ции, выполненная с использованием специализированных научных программных средств поддержки плани- 
рования эксперимента. 
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Введение. В научной и инженерной практике часто возникают задачи планирования экспери- 
мента. Чтобы эффективно спланировать, провести эксперимент, а также обработать его резуль- 
таты, необходимо воспользоваться специализированными программными средствами (ПС). Обыч- 
но применяются МАТЬАВ [1], ЗТАТТЬТТСА [2], Безюап-Ехрем [3] и др., но перечисленные про- 
граммные продукты имеют высокую стоимость, закрытый исходный код, а также не все являются 
кроссплатформенными, что затрудняет их применение на практике. Рассматриваются возможно- 
сти эффективной комбинации свободных научных библиотек Ру{Поп и созданного автором ПС ЕО- 
биррои [4] на примере решения задачи экспериментального поиска оптимальных параметров од- 
ного генетического алгоритма (ГА) аппроксимации. 

Описание исследуемого ГА. Рассматриваемый ГА предназначен для поиска наиболее длинного 
прямолинейного участка на заданной кривой и аппроксимации его отрезком (рис.1). 
В терминологии ГА каждый отрезок — особь, 
содержащая четыре гена: длина отрезка, угол 
наклона, координаты х и у центра отрезка. 


Скрещивание — обмен между двумя особями 
одним случайно выбранным геном. Мутация 
состоит в добавлении небольшой случайной 
величины к одному, произвольно выбранному, 
гену особи. 

Постановка задачи. Скорость поиска реше- —14 
ния является критичной характеристикой при 
практическом применении исследуемого ГА. На Рис.1. Пример результата работы ГА 

нее существенно влияют значения множества (ни обоего отека отронирио 
параметров алгоритма, важнейшими из кото- 

рых являются вероятности мутации и скрещивания особей. Разные значения этих двух парамет- 
ров (при фиксированных остальных) приводят к различной скорости сходимости ГА 
к искомому решению (табл.1). 
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Таблица 1 
Количество итераций ГА, необходимых для достижения решения, 
близкого к оптимальному, при различных значениях вероятностей мутации и скрещивания 




















№ п/п Вероятность мутации Вероятность скрещивания Кол-во итераций 
1 0,5 0,5 24 
2 0,85 0,3 28 
3 0,4 0,4 26 
4 0,6 0,4 25 
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В связи с этим возникает задача поиска оптимальных или близких к ним значений конфи- 
гурационных параметров ГА, обеспечивающих наибольшую скорость сходимости алгоритма. По- 
скольку ГА является разновидностью эвристических алгоритмов, решение данной задачи возмож- 
но только методом экспериментального исследования. 

Определение условий эксперимента. В качестве управляемых факторов х, и х› выбраны 


значения соответственно вероятностей мутации и скрещивания. Другие основные параметры ГА: 
размер популяции — 2000 особей; отношение размера селекции к размеру популяции - 0,6; ис- 
пользуемый алгоритм селекции - турнирный отбор с численностью турнира А =2. 

В качестве критерия оптимизации, который должен отражать скорость сходимости алго- 

ритма, принято количество итераций, необходимых для получения решения, близкого к опти- 
мальному. Количество итераций берется с отрицательным знаком для приведения к задаче мак- 
симизации функции отклика. 
Восхождение по поверхности отклика. Поставлен первый полный факторный эксперимент 
(ПФЭ) в области факторного пространства, обозначенной 1 на рис.2. В каждой точке плана про- 
водилось по 20 параллельных опытов. Полученная линейная модель у=-26,55 - 3,45х/ + 0,8х., 
предоставила информацию о направлении роста функции отклика. Выполнено несколько шагов 
восхождения по поверхности отклика (пунктирная линия на рис.2)) в сторону градиента функции. 
Вычисление коэффициентов модели, проверка ее адекватности, а также расчет траектории вос- 
хождения выполнялись средствами ПС ЕОбирром. 
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Рис. 2. Расположение областей экспериментов 
в факторном пространстве 


Для автоматизации эксперимента использовались скрипты, написанные на языке РУ{оп, с 
применением библиотеки для работы с матрицами МитРу [5], что позволило производить вычис- 
ления со скоростью, не уступающей скорости мощной вычислительной среды МАТЕАВ [6]. 

В области, обозначенной 2 на рис.2, был поставлен второй ПФЭ и получена новая линей- 
ная модель у = -—23,19 —0,058х, + 0,81х. . Снова были выполнены шаги восхождения по поверхно- 


сти отклика, в результате чего была достигнута область предполагаемого оптимума. 

Планы для всех экспериментов были сгенерированы при помощи ПС ЕОбуррок, которое 
предоставляет пользователю удобный интерфейс для выбора количества факторов, интервалов 
их варьирования и типа генерируемого плана эксперимента. На рис.3 представлено диалоговое 
окно выбора параметров плана. 
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Структура регрессионной модели 
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Рис.3. Окно программы ЕО5иррок для подготовки ПФЭ 


Описание экстремальной области. В силу того, что достигнутая подобласть факторного про- 
странства (обозначенная 3 на рис.2) практически вплотную приблизилась к допустимой границе, 
было решено провести заключительный эксперимент на основе центрального композиционного 
плана (ЦКП). Для повышения точности результата было увеличено количество параллельных 
опытов до 30. В табл.2 представлена использованная матрица планирования, а также результаты 


ОПЫТОВ (где у - среднее значение отклика по параллельным опытам). 
































Таблица 2 
Матрица планирования ЦКП и результаты проведенных измерений 

[е] == [е] — 
№ хм Хх У ыы Хх х› У 

1 1 1 23,73 9 0 0 22,57 
2 —1 1 21,87 10 0 0 21,2 
3 1 —1 23,53 11 0 0 23,03 
4 —1 —1 22,27 12 0 0 23,17 
5 —1,369 0 -27,07 13 0 0 21,33 
6 1,369 0 -24,73 14 0 0 21,7 
7 0 —1,369 23,63 15 0 0 22,33 
8 0 1,369 21,77 = № — = 


























Анализ результатов. Полученная в результате эксперимента ЦКП квадратичная модель после 
удаления незначимых коэффициентов в кодированных величинах факторов имеет следующий 


вид: 


у=-22,13 + 0,36х, — 0/15х,х, —1,54х,*. 
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График данной модели представлен на рис.4. Все графики данной статьи были построены 
при помощи специализированной библиотеки Маро [7] для языка Руфоп, имеющей богатые 
возможности в области визуализации данных. 
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Рис.4. График квадратичной модели поверхности отклика 
в области в натуральных величинах факторов 


Из анализа полученной квадратичной модели следует, что оптимальные значения вероят- 
ностей мутации и скрещивания находятся в следующих интервалах: х„: 0,35-0,4; х,: 0,9-1,0. 
Серия дополнительных опытов показала, что значение функции отклика в данной области соот- 
ветствует в среднем 22 итерациям, что является объективным улучшением результатов (табл.1). 
Выводы. Были исследованы свойства специализированного ГА аппроксимации, определены 
близкие к оптимальным значения его основных параметров, позволяющие получить наибольшую 
скорость сходимости. Показана эффективность применения в научной деятельности свободных 
программных средств на основе языка Ру{Поп в сочетании с ПС ЕОбуррок. 
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